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Abstract 

We obtain the horizontal air shower rate from the (EeV) high energy 
neutrino flux predicted in some topological defect scenarios as the source for 
the highest energy cosmic rays. Emphasis is made on the different character 
of the events depending on the neutrino flavor an type of interaction. We 
show that the bound for muon poor showers in the 10 5 — 10 7 GeV energy 
range is violated by maximal predictions for superconducting cosmic string 
neutrino fluxes, we compare it to other neutrino flux limits and we discuss 
the future of such measurements to further constrain these models. 
The establishment of cosmic rays of energies above 10 20 eV [|l| and the prospects for 
the construction of a giant air shower array for their study, the Pierre Auger project ||, 
has stimulated great activity in the search of theoretical models for the production of these 
particles. The annihilation of topological defects (TD's) has been recently proposed as 
the origin of the highest energy cosmic rays || although there is some controversy about 
this possibility ||]. This mechanism avoids the main difficulties in the conventional shock 
acceleration and would imply large fluxes of gamma rays and neutrinos up to energies in 
the EeV range and well above ||. Active Galactic Nuclei (AGN), the most powerful known 
objects, have also been suggested as the possible origin of the highest energy cosmic rays ||. 
If the very high energy gamma rays detected from AGN @ come from pion decay, in proton 
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acceleration models, they must also emit high energy neutrinos. Models for acceleration in 
the AGN jets predict neutrinos also extending to the EeV range PHTU|. 

Some of these models have already been shown to predict secondary photon and neu- 
trino fluxes in conflict with experiment. Gamma ray searches (~ 100 MeV) constrain certain 
combinations of extragalactic magnetic fields and TD models because of electron positron 
pair production in the magnetic field [JO]. Also the Frejus underground muon detector has 
not observed the muon neutrino prediction for Superconducting Cosmic Strings (SCS) jl2| . 
Clearly there are many uncertainties in these models such as the normalization, the mass 
of the X particle in the GUT scale, Mx, the extragalactic magnetic field, the fragmentation 
functions involved as well as the propagation of each of the produced particles in the extra- 
galactic magnetic field and background radiations, taking into account the evolution of the 



Universe. It has been claimed by different groups ||13l , |T4|| that tuning the free parameters 
in the model (Mx, the maximum of the injection spectrum and the extragalactic magnetic 
field), it is still possible to explain the UHE cosmic events, avoiding the 100 MeV gamma 
ray constraints ||T^| . 

Large scale high energy neutrino detectors are expected to further test these models || . 
As the earth becomes opaque for neutrinos already in the 100 TeV to 1 PeV range, EeV 
neutrinos can only be detected for zenith angles between vertical down going and horizon- 
tal. In such case the experiment must have some control over the energy of the muon to 
separate the neutrino signal from that of atmospheric muons which poses a serious back- 
ground. Underground muon experiments rely on their overburden for this purpose. For high 
energies they must search for close to horizontal neutrinos where the overburden and energy 
uncertainties are largest. Detectors in water or ice may, in addition, detect Cerenkov light 
from neutrino showers with the advantage of being also sensitive to the electron neutrino. 
They require good directional and shower reconstruction capabilities to separate neutrino 
showers from those induced by muons. Horizontal Air Shower (HAS) measurements provide 
a third alternative as they are due to high energy penetrating particles such as muons and 



neutrinos |16| , including electron neutrinos, and shower size determinations can be used as 



an energy threshold. Existing air shower arrays can constrain neutrino fluxes [17,18| and 



their relevance for direct neutrino detection will become more important when the Pierre 



Auger project is constructed [19,20 



In this paper we calculate the HAS rate due to neutrino fluxes from TD and recent 
AGN models separating explicitly the muon poor HAS in order to compare it to the bound 
published by the AKENO group [^TJ. The v e + v e flux prediction in the maximal super- 
conducting string model violates this bound. We compare the different shower channels in 
this model showing that about half the total v e + v e contribution is muon poor. The muons 
produced by charged current muon neutrino interactions in the atmosphere dominate the 
"standard" atmospheric muon flux at high energies and we calculate the HAS contribution 
from secondary bremsstrahlung from these muons. This is the dominant contribution to 
muon poor showers from + v^. We finally compare the AKENO bound to other bounds 
emphasizing the differences in energy range and we briefly address the potential capabili- 
ties of typical existing and planned air shower arrays to detect neutrino induced horizontal 
showers. 

Horizontal Air Showers. Extensive air showers initiated by primary cosmic rays are 
strongly suppressed at large zenith angles because the slant depth of the atmosphere rises 
from ~ 30 radiation lengths in the vertical direction to ~ 1,000 in the horizontal direction 
and the absorption depth of showers only rises logarithmically with primary energy. It has 
been known for long that only weakly interacting particles, such as muons and neutrinos, 
may penetrate such large depths to interact and produce a shower sufficiently close to be 
detected at ground level. HAS have been observed in the 60's by several groups for shower 
sizes between 10 3 -10 5 particles p2| . The shower rate observed is consistent with conventional 



muon production by hard bremsstrahlung [^,^,|T^,Q . For shower sizes above 10 5 particles, 
AKENO has published an upper bound on muon poor air showers at zenith angles greater 
than 60° [^TJ. These results have been used as a bound for the high energy muon fluxes 
induced by the semileptonic decay of charm particles produced by primary cosmic rays in 
the atmosphere [|17], and as a non-accelerator hint on the very high energy charm cross 



section |25| . 



There are multiple shower channels for neutrino showers. Deep Inelastic Scattering (DIS) 
interactions of the neutrino with air nucleons via charged and neutral currents produce a 
hadronic shower of energy yE u , where y is the energy fraction transferred to the nucleon. In 
charged current electron-neutrino (antineutrino) interactions the emitted electron (positron) 
generates an electromagnetic cascade carrying the remaining energy of the initial neutrino 
(1 — y)E v that is superimposed to the cascade initiated by the hadron debris. For muon 
neutrinos the charged current interaction produces a secondary flux of high energy muons 
which could in turn induce horizontal showers further into the atmosphere in the same way 
the atmospheric muons do. The interaction of neutrinos with atomic electrons is in general 
suppressed by the ratio of the electron to the proton rest mass except for the Glashow 
resonance, V e + e~ — *> W~, which dominates over all other processes at the resonance energy 
Ep e ~ 6.4 10 6 GeV p6 |. The decay of the W boson into an electron and an antineutrino 



generates a purely electromagnetic shower with energy E ~ 3 10 6 GeV. Its decay in the 
muon channel will also contribute to the horizontal muon flux at these energies, while the 
tau channel will in turn decay into hadrons 64% of the times and into a muon or an electron 
18% of the times each. 

Following ref. [|17] the differential rate of horizontal cascades at fixed shower size, N e , 
is obtained similarly for neutrinos and muons integrating the parent flux, the y differential 
cross section of the corresponding interaction and the atmospheric depth: 

MN e ) = J o dtJ Q ^ = —— j __ , (1) 

Here y is y, 1 — y or 1 depending on the interaction that produces the shower. For muon 
bremsstrahlung the energy of the photon which originates a shower at a given depth and of 
a given shower size is determined, on average, by the inverse of the Greisen parametrization 
E(N e , t) P?fl . The energy of the parent muon is given by E/y. The factor dE/dN e is just the 
Jacobian of the transformation. For neutral current interactions the energy of the hadron 



debris is inferred inverting the parametrization of hadronic showers due to Gaisser [p8[| . In 



the case of charged current electron neutrino interactions, eq. (|T|) applies with y = 1, because 
all the energy is transferred to the shower; and a more complicated shower size depending 
on y and combining both parametrizations has been used. 

We have performed the previous calculations using the neutrino flux predictions from 
topological defects [^j and from AGN models ||, and we have assumed that the fluxes for 
u e , z/ e , and z/ M are in the ratio 1:1:2:2. In fig.(0) we show the + differential fluxes 
in these models compared to atmospheric neutrinos. For the neutrino cross section we have 
used the MRS(G) parton distribution functions f3(| extrapolated to low x as described in 
|19| . We have checked that our results are not sensitive to using alternative parametrizations 
for shower size such as that of Fenyves et al. [pi , nor to the multiplicity of hadrons in the 



DIS interactions. In fig.(|2]) we show the total integral horizontal air shower rates for the 
TD and AGN models as a function of the shower size, compared to Tokyo data and to 
expectations from atmospheric muons. 

Muon poor showers. The AKENO bound on horizontal air showers has been obtained 
on the basis of searching for muon poor showers. This has allowed to search for showers 
of zenith angle already above 60°; but it means that only the muon poor fraction of the 
HAS produced by a hypothetical neutrino flux has to be compared with this bound. Only 
the resonant channel decaying into an electron and an electron-antineutrino produces a 
purely electromagnetic shower directly in the interaction. Electron-neutrino charged current 
interactions in which the hadronic component carries less than the 20% of the initial energy 
E u (y < 0.2) produce muon poor showers according to AKENO requirements. Since the 



average < y > at high energies approaches 0.2 [[26[] a good fraction of the electron neutrino 
showers are muon poor. Both these contributions are calculated as described in the previous 
section. For the SCS model we compare the total rate of HAS, to that from electron neutrinos 
and antineutrinos only, separating the DIS and the resonant contributions in fig. ([$]). We 
also plot in the same figure the results for muon poor horizontal showers as described above 
separating the resonant and DIS channels as well. Of all the electron neutrino HAS, about 
50% are muon poor at high energies. 



Any channel for producing muons can also give pure electromagnetic showers through 
secondary muon bremsstrahlung or electron positron production. In spite of the suppression 
due to a double process its contribution can be relevant for muon poor showers at large zenith 
angles. We consider muons produced in charged current muon-neutrino interactions and 
those produced in the resonant interaction decaying into a muon and a muon antineutrino, 
these muons carry most of the energy of the interaction. Muons produced in hadronic showers 
or tau decays have lower energy and can be neglected in comparison to other uncertainties. 

For incoming muon neutrinos, we have firstly calculated the muon flux produced by the 
charged current interaction of horizontal muon neutrino fluxes through the 36000 g cm~ 2 
slant depth. This flux will give rise to purely electromagnetic horizontal showers due to hard 
processes of energy loss such as bremsstrahlung and electron positron production. These 
showers will generally develop well after the first shower initiated by the nucleon debris is 
absorbed. The calculation of this channel uses eq. ([1]) with the secondary muon flux, the 
muon cross section and y = y. We consider the dominant contribution, bremsstrahlung, 
and we use the cross section given by |3"2" |. The shower rate calculated in this way at 90° 



represents approximately 50% of the "direct" electron neutrino showers with a cut in y < 0.2. 
Secondary bremsstrahlung is the most important channel for muon poor showers induced 
by muon neutrinos and has a strong dependence on zenith angle. The result of secondary 
muon bremsstrahlung HAS in the SCS model, averaged over the zenith angles from 60° to 
90°, is also compared to the other channels in fig.([3|). However the contribution represents 
less than a 10% effect on the total number of muon poor showers because of this averaging. 

In fig. (H) we display the results for muon poor showers in the considered models together 
with the AKENO bound where it is shown how the model with p = 0, corresponding to 
superconducting cosmic strings, is ruled out by an ample margin. The figure also illustrates 
how the model with different time evolution (p = 0.5) is not far from being tested with 
current experiments assuming improvements in statistics. Clearly searching for HAS of any 
type should increase the expected signal, but care has to be taken to insure HAS are not due 
to ordinary cosmic rays. We have include in fig.(@) a rough theoretical estimate of the 90% 
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confidence upper limit for no observations during a year, assuming no ordinary cosmic ray 
showers penetrate at angles above 75°, for ideal air shower arrays of three areas ~ 10 5 m 2 , 
~ 10 2 km 2 and ~ 10 4 km 2 . These are typical order of magnitude sizes of existing and 
planned arrays. 

Limits. The information provided by Frejus, HAS and Fly's Eye, ruling out the same 
types of models, is complementary in the sense that they relate to very different parts of the 
neutrino spectrum and because HAS apply to electron neutrinos as well. It is interesting to 
compare the limit derived from the AKENO bound with other limits for high energy neutrino 
fluxes. At lower energies (~ 3 TeV), the Frejus limit |E| rules out the TD prediction with 



p = 0. The Fly's Eye limit ||18|| , also rules out this prediction at much higher energies 
~ 10 5 TeV. The limit from AKENO lies in between, in the PeV region. The EAS-TOP 



collaboration |33] has also published a less constraining limit on high energy neutrino fluxes 
from horizontal showers also near the PeV region. In fig. ([3]) we show the flux predictions 
together with these limits. This figure, however, must be interpreted with care as these 
limits can be very flux- dependent. We illustrate this by plotting the limits corresponding to 
constant spectral index fluxes with different indices (7 = 1.1,2,3). In this respect it would 
be useful that all experiments measuring Air Showers such as EAS-TOP, CYGNUS and 
Fly's Eye (measuring upcoming showers), gave their results as air shower rates like AKENO 
as this would simplify the comparisons between experiments and would be a step forward 
to the establishment of these results. 

To further illustrate this point we have considered a simple neutrino flux with constant 
spectral index (7) from the GeV to the EeV region <p(E) = AE~ y . For each of the above 
limits there is a curve in the 7— A parameter space plane. This is shown in Fig. There are 
several interesting things to point out from this figure. The best limit corresponds to Fly's 
Eye for small spectral indices (as is the case for the TD models), and to deep underground 
experiments such as Frejus if the index is large. Besides, for 7 ~ 2 all experiments have 
similar performances within an order of magnitude. 

HAS (within which we include Fly's Eye results) are a tool for constraining TD scenarios 



via the neutrino-induced showers produced along with the UHE cosmic rays and they rule 
out superconducting cosmic string models normalized to the highest energy cosmic ray 
spectrum. The prospects of running the largest existing air shower arrays using horizontal 
triggers, in order to get not only the electromagnetic shower rate bounds, but the rate of 
all Horizontal Air Showers, should further constrain the present TD models. Clearly the 
possibility of searching for HAS with future air shower arrays such as the Pierre Auger 
Project, with an area about 60 times larger than AKENO, if viable, would undoubtfully 
provide a much stronger constraint on these models. 

Acknowledgements 

We thank Jaime Alvarez Muhiz, Francis Halzen, Karl Mannheim, Gonzalo Parente and 
Giinter Sigl for helpful discussions. This work was partially supported by CICYT under 
contract AEN96-1773, by the Xunta de Galicia under contract XUGA 20604A96 and by the 
EC (R.A. V.), under contract ERBCHBICT941658. 



8 



REFERENCES 

[1] D.J. Bird et al, Phys. Rev. Lett. 71 (1993) 3401; M. Nagano et al., J. Phys. G: Nucl. 
Part. Phys. 18 (1992) 423. 



[2] 
[3] 



[4] 
[5] 
[6] 
[7] 
[8 
[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 



The Pierre Auger Project, Design Report, FERMILAB Report, Oct. 1995. 

J. P. Ostriker, C. Thompson and E. Witten, Phys. Lett. B180 (1986) 231; C.T. Hill, 
D. N. Schramm, and T.P Walker, Phys. Rev. D36 (1987) 1007; P. Bhattacharjee and 
N.C. Rana, Phys.Lett. B246 (1990) 365. 

A.J. Gill and T.W.B. Kibble, Phys. Rev. D50 (1994) 3660. 

G. Sigl, D.N. Schramm and P. Bhattacharjee, Astropart. Phys. 2 (1994) 401. 

For a review, see F. Halzen, T.Stanev, T.K. Gaisser, Phys. Rept. 258 (1995) 173. 

D.J. Thompson et al, Ap.J. Suppl. 101 (1995) 259. 

K. Mannheim, Astropart. Phys. 3 (1995) 295 

K. Mannheim, in: Workshop on High Energy Neutrino Astrophysics, Aspen (1996). 
R.J. Protheroe; Ibid. 

R.J. Protheroe and T. Stanev, preprint |astr-ph / 9605036 . 



W. Rhode, K. Daum, et al. Proc. 24th ICRC, Rome 1995. vol. 1, p. 781. Also, W. 
Rhode, et al. Astropart. Phys. 4 (1996) 217. 



R.J. Protheroe and P.A. Johnson, |astro-ph/9605006| . See also R.J. Protheroe and P.A 



Johnson, Astropart. Phys. 4 (1996) 253. 
G. Sigl, S. Lee and P. Coppi, Fermilab-Pub-96/087-A, |astro-ph/9604093 



X. Chi et al. Astropart. Phys. 1 (1993) 239; X. Chi et al. Astropart. Phys. 1 (1992) 129. 
V.S. Berezinsky and A. Yu. Smirnov, Astrophys. Space Science 32 (1975) 461. 



[17] E. Zas, F. Halzen and R. A. Vazquez, Astropart. Phys. 1 (1993) 297. 
[18] R.M. Baltrusaitis et al. Phys. Rev. D31 (1985) 2192. 

[19] G. Parente and E. Zas; in Proceed, of the 7th Int. Symp. on Neutrino Telescopes, pp 
345, Ed. M. Baldo Ceolin, Venice (1996). 

[20] J. Capelle, J.W. Cronin, G. Parente and E. Zas; in preparation (1996). 

[21] M. Nagano et al, J. Phys. G: Nucl. Part. Phys. 12 (1986) 69. 

[22] M. Nagano et al, J. Phys. Soc. Japan. 30 (1971) 33; S. Mikamo et al, Lett, al Nuov. 
Cim. 34 (1982) 273. 

[23] P. Kiraly, M.G. Thompson and A.W. Wonfendale, J. Phys. A: Gen. Phys. 4 (1971) 367. 

[24] G. Parente, A. Shoup and G.B. Yodh, Astropart. Phys. 3 (1995) 17. 

[25] M.C. Gonzalez-Garcia, F. Halzen, R.A. Vazquez and E. Zas, Phys. Rev. D49 (1994) 
2310. 

[26] R. Gandhi et al, Astropart. Phys. 5 (1996) 81. 

[27] K. Greisen, Prog, in Cosmic Ray Phys., ed. J.G.Wilson, vol. Ill p. 1, North Holland 
Publ. Co. (1993). 

[28] T.K. Gaisser, Cosmic rays and Particle physics, Cambridge University Press (1990) 279 
P- 

[29] P. Bhattacharjee, C.T. Hill, D.N. Schramm, Phys. Rev. Lett. 69 (1992) 567. 
[30] A.D. Martin, W.J. Stirling and R.G. Roberts, Phys. Lett. B354 (1995) 155. 
[31] E.J. Fenyves et al, Phys. Rev. D37 (1988) 649. 

[32] S.R. Kelner, R.P Kokoulin and A. A. Petrukhin, Moscow preprint 024-95, 1995. 
[33] M. Aglietta et al. EAS-TOP coll. Proc. 24th ICRC, Rome 1995, Vol. 1, p. 638. 

10 



FIGURES 



11 



10° 



EAS-TOP 




E„ (GeV) 



12 



FIG. 1. Differential (f M + v^) fluxes of ref. 29 for p = (solid line) and for p = 0.5 (dotted). 
The dashed line is the prediction of ref. 9 and the dot-dashed line is the atmospheric flux. Also 
shown are some published limits as marked. The three parallel lines illustrate the change in the 
limits for different spectral indices (1.1, 2 and 3 from bottom to top). 



13 



,-13 



Tokyo data 
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FIG. 2. Integral horizontal shower rate as a function of shower size (Ne) for different models. 
TD neutrinos with p = (solid line), p = 0.5 (dotted), AGN neutrinos (dashed) and atmospheric 
muons (dot-dashed). Diamonds represent the HAS data from the University of Tokyo. Horizontal 
lines marked as A, B and C represent the expected limit of non observation in a year made with 
ideal detectors of areas 10 5 m 2 (A), 10 8 m 2 (B), and 10 10 m 2 (C). See text. 
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FIG. 3. Partial HAS channels for the TD model with p = 0. From top to bottom: total rate, 
{y e +v e ) DIS contribution, the muon poor [v e -\-v e ) DIS contribution, the total resonant contribution 
ipe), the secondary muon bremsstrahlung and the resonant contribution into the electron channel. 
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FIG. 4. Muon poor integral shower rate for the three models. TD p = (solid line), p = 0.5 
(dotted) and AGN (dashed). 
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FIG. 5. Limits on the A — 7 space (see text) given from different experiments, as marked in 
the figure. The energy attached to each line gives an indication of the energy involved in the 
corresponding experiment. The dotted line (marked as KGF) corresponds to the point source limit 
(Markarian 421) from the Kolar Gold Field experiment. 
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